arXiv:1505.05436vl [astro-ph.IM] 20 May 2015 


Mon. Not. R. Astron. Soc. 000,[H[j0](2014) 


Printed 21 May 2015 (MN style file v2.2) 


Optical Turbulence Characterization at LAMOST Site: 
Observations and Models 


L.-Y. Liu 1 ’ 2 *, C.Giordano 2 , Y.-Q. Yao 1 , J.Vernin 2 , M. Chadid 2 , H.-S. Wang 1 , 

J. Yin 1 , Y.-P. Wang 1 

1 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China 

-Universite de Nice-Sophia Antipolis, Observatoire de la Cote d’Azur, CNRS-UMR7293, Laboratoire Lagrange, 06108 Nice Cedex 2, France 


Accepted. Received; in original form 


ABSTRACT 

Atmospheric optical turbulence seriously limits the performance of high angular resolution 
instruments. An 8-night campaign of measurements was carried out at the LAMOST site in 
2011, to characterize the optical turbulence. Two instruments were set up during the cam¬ 
paign: a Differential Image Motion Monitor (DIMM) used to measure the total atmospheric 
seeing, and a Single Star Scidar (SSS) to measure the vertical profiles of the turbulence Cjdi) 
and the horizontal wind velocity V{h). The optical turbulence parameters are also calculated 
with the Weather Research and Forecasting (WRF) model coupled with the Trinquet-Vernin 
model, which describes optical effects of atmospheric turbulence by using the local meteoro¬ 
logical parameters. This paper presents assessment of the optical parameters involved in high 
angular resolution astronomy. Its includes seeing, isoplanatic angle, coherence time, coher¬ 
ence etendue, vertical profiles of optical turbulence intensity Cjih) and horizontal wind speed 
V{h). The median seeing is respectively 1.01 arcsec, 1.17 arcsec and 1.07arcsec as measured 
with the DIMM, the SSS and predicted with WRF model. The history of seeing measure¬ 
ments at the LAMOST site are reviewed, and the turbulence measurements in this campaign 
are compared with other astronomical observatories in the world. 

Key words: Turbulence - Atmospheric effects - Instrumentation: high angular resolution - 
Site testing 


1 INTRODUCTION 

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope 
(LAMOST), is a special quasi-meridian reflecting Schmidt tele¬ 
scope; its particular design allows both a large aperture (effective 
aperture of 3.6 - 4.9 meter) and a wide field of view (FOV=5°) 
(Richard 2008, Cui et al. 2012). The main scientific goal of LAM¬ 
OST focuses on the spectroscopic survey in the Galaxy and ex- 
tragalactic (Zhao et al. 2012). LAMOST located at Xinglong ob¬ 
servatory of National Astronomical Observatory of China(NAOC), 
which has been regarded as one of the important optical astronomi¬ 
cal observatories in China. The Xinglong Observatory is located in 
Hebei province, 170 kilometers northeast of Beijing, where many 
telescopes are operating including the LAMOST, the 2.16 meter 
optical telescope and several 1 meter telescopes. A site testing cam¬ 
paign was performed in 2011 to obtain the atmospheric optical pa¬ 
rameters at Xinglong observatory, for the application of Adaptive 
Optics(AO), the active optics and the high resolution spectrometer 
on the LAMOST and the 2.16 meter telescope(Zhang et al. 2004, 
Su et al. 2012). 
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The early work on optical turbulence characterization at Xin¬ 
glong observatory was carried out by Wu et al. (1996) and Song 
et al. (1998). They carried out a 4 night campaign of multi¬ 
instrument measurements in December 1994. The integrated see¬ 
ing was measured with a DIMM, the surface layer turbulence 
was detected with microthermal sensors and an acoustic radar. 
Based on the photometry dataset in 1995-2001 of Beijing-Arizona- 
Taipei-Connecticut(BATC) telescope, Liu et al. (2003) analyzed 
the Polaris images by measuring the Full Width at Half Maxi- 
mum(FWHM) to understand seasonal features of seeing condition. 
Liu et al. (2010) performed DIMM seeing measurements on the 
top of LAMOST building, in order to optimize the slit width of 
LAMOST spectrometer. Yao et al.(2012) analyzed the seeing con¬ 
ditions with the BATC photometric data in 1995-2011, and their 
results could be referred as the long-term (more than 15 years) see¬ 
ing variation at Xinglong observatory. 

In this paper we present the results of the site testing campaign 
at the LAMOST site in 2011. In this campaign, two kinds of instru¬ 
ments were employed in order to obtain turbulence profiles and in¬ 
tegrated atmospheric parameters simultaneously. We have also cal¬ 
culated the atmospheric parameters at the same time with a numer¬ 
ical method based on the WRF model coupled with the Trinquet- 
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Vernin model (Giordano 2013, Trinquet 2007). In this paper we do 
not intend to show the reliability of WRF model to foresee a 3D 
map of the optical turbulence due to too few days of observation. 
This goal was already demonstrated by Giordano (2013). 

In Section [2] we briefly recall the optical turbulence param¬ 
eters referring to high angular resolution techniques, and describe 
the DIMM, SSS, and WRF model used in this campaign. Section[3] 
presents the results obtained in the campaign. The comparison and 
discussion of the results are given in Section[4] in which we review 
the measurements at LAMOST site, and compare with other sites 
in the world. Section[5]is a summary. 


2 INSTRUMENTS AND MODELS 

Atmospheric optical turbulence is the main factor of the degrada¬ 
tion of image resolution. The major characteristics of optical turbu¬ 
lence are described as profiles and integrated parameters. The two 
main vertical profiles are the profile of the refractive index structure 
constant C? t (h), and the profile of the horizontal wind velocity V(h). 

In order to improve the techniques of high angular resolution 
(AO, spectroscopy, etc.), one needs knowledge of the atmospheric 
integrated parameters, such as, the Fried’s radius ro or the seeing so, 
the coherence time r 0 , and the isoplanatic angle 6q. The integrated 
parameters can be retrieved by the profiles above mentioned. 

The Fried's radius is: 


ro 


0.185T 675 



r 3/5 

ClQi)dh\ 


( 1 ) 


where ,1 is the wavelength, h is the altitude above the ground. 

The seeing is one of the most important turbulence parameters 
to astronomy, and usually the seeing value of observatory sites is 
around 1 arcsec in the visible band. As usual, here we assume ,1 = 
0.5yum. 

Seeing e 0 (eq. 0. coherence time r 0 (eq. [3} and isoplanatic 
angle do (eq.0, can be derived from the vertical profiles of Cf t (h) 
and wind speed V(h), as follows: 


so = 0.98T/ro, 


( 2 ) 



Figure 1. The DIMM and SSS instruments are installed in the Public Dome, 
Xinglong observatory. 


Assuming a Kolmogorov model for the optical turbulence, the see¬ 
ing can be calculated from the variance of the differential motion 
between two star images. Therefore, DIMM seeing monitor is al¬ 
most insensitive to tracking errors, and the instrument yields the 
integrated seeing from telescope pupil to the top of the atmosphere. 
Nowdays, the DIMM technique is universally employed in site test- 
ingfTokovinin, 2002; Aristidi, 2005). 

According to the principle mentioned above, our DIMM in¬ 
strument has the following configurations: 

(i) a Meade LX200GPS telescope, with the entrance pupil of 20 
cm, and the focal length of 200 cm; 

(ii) a mask with two sub-aperture, with diameter of 5 cm, and a 
separation of 15 cm; 

(iii) an optical wedge installed on one sub-aperture, to adjust an 
appropriate separation between the two star images at focus; 

(iv) a Lumenera SKYnyx 2.0M CCD camera, attached to 
the telescope focal plane for fast sampling star images (At = 
0.5/ni),with a 640 x 480 format and the pixel size of 7.4 gm; 

(v) two computers, one calculating seeing values and closed- 
loop guiding, and the other remotely monitoring the status via in¬ 
ternet. 


T 0 = 0.058T 675 


do = 0.058T 675 



r 3/s 

\V(h)\ 5/3 C 2 „(h)dh\ 


h 5n C*(h)dh 


- 3/5 


(3) 

(4) 


The coherence etendue Go (Lloyd et al. 2004), is a comprehen¬ 
sive evaluation for adaptive optics with the combination of Fried’s 
radius r 0 , isoplanatic angle d 0 and coherence time r 0 , defined by the 
following formula: 


G 0 = TqTo^q. 


(5) 


2.1 DIMM 

The DIMM is a small transportable instrument, with the differential 
technique to accurately measure the seeing conditions, as described 
by Sarazin(1990) and Vernin(1995). A small telescope is equipped 
with a mask made of two sub-apertures, in order to obtain simulta¬ 
neously two separated images of the same star on the focal plane. 


2.2 Single Star Scidar 

The Single Star Scidar is a new member in Scidar family (Vernin 
& Roddier 1973; Fuchs, Tallon & Vernin 1998; Avila et al. 1998; 
Habib et al. 2006; Garcfa-Lorenzo & Fuensalida 2011; Aviles et al. 
2012; Masciadri et al. 2012; Shepherd et al. 2014). This instrument 
can deliver the distribution of C\ with height, by using a small tele¬ 
scope. The SSS technique analyzes the spatial-temporal auto- and 
cross-correlation functions of less than 1 ms exposure images from 
a single star, which allow to retrieve the vertical profiles of opti¬ 
cal turbulence intensity C\{\i) and horizontal wind speed \(h) from 
the ground to the top of the atmosphere, yielding the determina¬ 
tion of seeing, isoplanatic angle and coherence time (See Section 
0. The SSS technique in detail can be referred to Habib(2006) and 
Vernin(2009). In order to obtain the contribution of surface layer 
turbulence, the SSS also employs the so-called ’’Generalized Sci¬ 
dar” technique (Fuchs et al. 1998). Although the SSS is a low- 
resolution optical turbulence profiler, it needs only 40 cm aperture 
telescope, showing a practical advantage of being transportable for 
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field site testing (Vernin et al. 2011, Giordano et al. 2012, Liu et al. 
2012, Liu et al. 2013a). 

In order to obtain optical turbulence profile at the LAMOST 
site, we have constructed a new SSS with commercial components 
(Liu et al. 2013b). The main specifications of our SSS are the fol¬ 
lowing. 

(i) a 40 cm Meade M16 tube, on an Astro-Physics 1200 equato¬ 
rial mount, with the focal length of 400 cm; 

(ii) a collimating lens after the focus of the telescope, with 10 
mm focal length, to make the beam parallel; 

(iii) a CCD camera, after the collimating lens for fast sampling 
the star scintillation pattern. The fast readout Pixelfly CCD-200XS, 
640x480 format with pixel size of 9.9x9.9 /rnr, allows continuous 
acquisition with a high frequency rate without loss of any image. 
The exposure time is usually taken as 1 ms, every 5.6 ms, the pixel 
scale is 0.51 arcsec/pixel; 

(iv) an auto-guiding system, made with a cube splitter; 

(v) two control computers, one capturing scintillation images 
from the CCD camera and computing the spatial auto-correlation 
and cross-correlation images, the other serving as real time track¬ 
ing through a CCD guider. 

2.3 VVRF Model 

Some models have been developed in order to estimate the op¬ 
tical turbulence parameters from meteorological and topography 
data (Coulman et al. 1988; Dewan et al. 1993; Masciadri et al. 
1999a; Masciadri et al. 1999b; Trinquet & Vernin 2007; Masci¬ 
adri et al. 2013). The Weather Research and Forecasting (WRF) 
model is a mesoscale non-hydrostatic numerical weather predic¬ 
tion system, used for both operational forecasting and atmospheric 
research. The WRF model allows to forecast the atmospheric pa¬ 
rameters at each point of a given 3D grid, such as the pressure P, 
the temperature T, the wind velocity components (u,v,w) and the 
relative humidity (Skamarock et al. 2008). 

The atmospheric optical parameters can be computed by the 
WRF model coupled with the Trinquet-Vernin model (Trinquet H. 
& Vernin J. 2006; Trinquet & Vernin 2007; Giordano et al. 2014). 
The Trinquet-Vernin model is a statistical model deduced from 160 
meteorological balloons analysis, and allows to compute the verti¬ 
cal profile of the Cj, from the vertical profiles of potential tempera¬ 
ture and wind velocity. A first study about the capability of WRF to 
predict the optical conditions is available in Giordano et al. (2013). 

For the optical turbulence characterization during the cam¬ 
paign at the LAMOST site, the model was initialized as the fol¬ 
lowing: 

(i) horizontal coarse grid, composed of lOOx 100 grid points,and 
spaced in Ax = Ay = 27 km; 

(ii) three consecutive nests with a ratio of 1/3, and therefore, the 
finest horizontal grid is Ax = Ay = 1 km; 

(iii) 87 vertical levels from 0km to 20km, with a higher resolu¬ 
tion within the low atmosphere A/?, = 50 m, and A/igy = 500 m; 

(iv) the NCEP FNL (Final) Operational Global Analysis data 
(http://rda.ucar.edu/! was used to initialize the grid, and to set the 
boundary conditions. 

The main parameters used for the simulation are: 

(i) the microphysics scheme used is the Thomson scheme 
(mp.physics = 8), 

(ii) the Rapid Radiative Transfer Model (RRTM) scheme is used 
for the long wave radiation (ra_lw_physics=l), 



Figure 2. The location of LAMOST on the Xinglong observatory (the photo 
was taken by Chen Yingwei). The DIMM and SSS instruments are located 
200 meters south-east of the LAMOST. The 2.16 meter telescope is located 
at the north-west of the LAMOST. 


Statistical distribution of e 0 with DIMM in 2011 



e Q [arcsec] 

Figure 3. The histogram (stairs) and cumulative histogram (dashed line) of 
DIMM seeing values at the LAMOST site by the site testing campaign in 
2011 . 

(iii) the Goddard shortwave scheme is used for the shortwave 
radiation (ra_sw_physics=2), 

(iv) the Mellor-Yamada-Janjic scheme is used for the planetary 
boundary layer (PBL) (bl_pbl_physics=2). 


3 SITE TESTING CAMPAIGN AND RESULTS 

The site testing campaign was carried out over 8 nights from April 
18 to September 26, 2011. The DIMM and SSS were installed 
in the Public Dome of Xinglong observatory, 4m above ground 
level, which is located about 200 meters south-east of the LAM¬ 
OST building, as shown in Fig|2] The Public Dome can be com¬ 
pletely open during observations. The two instruments, separated 
by 3 meters, were operated with the same target simultaneously. 
This configuration allows us to collect a substantial database and 
make cross-comparisons of the results. 

3.1 DIMM results 

The DIMM seeing measurements were carried out from April 18 
to 19 and from September 22 to 26. A total of 1738 seeing mea- 
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Figure 4. Histogram (solid line) and cumulative distribution (dashed line) of atmospheric parameters computed with the SSS measurements at the LAMOST 
site. The parameter distributions are seeing eq (top left), the isoplanatic angle do (top right), the coherence time To (bottom left), and coherence etendue Go 
(bottom right). 


surements has been analyzed. One DIMM seeing value was calcu¬ 
lated using a set of 100 images. The final seeing results have been 
corrected for the zenith angle. The observing targets for the DIMM 
seeing monitoring were oLeo (mV= 1.35 mag), a Boo (mV= -0.04 
mag), a Lyr (mV= 0.03 mag), a Aql (mV= 0.77 mag), /? Ori (mV= 
0.12 mag), a Cas(2.24 mag) and a Cyg (mV* 1.2 mag, variable 
star). 

Fig[3] shows the histogram and the cumulative distribution of 
the entire DIMM data set. The mean seeing value is 1.09 arcsec, and 
the median value 1.01 arcsec. There are 25% of the seeing values 
better than 0.9 arcsec, and 75% better than 1.2 arcsec. 

3.2 Single Star Scidar results 

During the campaign in 2011, a total of 6011 Cj,(h) and |v(/7)| pro¬ 
files were obtained with the SSS from April 18 to 19, May 13, and 
from September 23 to 26. The SSS spatiotemporal cross-correlation 
functions were obtained every 11.2 sec corresponding to a set of 
2000 images. Due to measurement errors be caused by continu¬ 
ously changing sky background, especially in the night sky with 
moon, the sky background changing with time was taken into ac¬ 
count, and checked every hour during the SSS measurement. The 


observing targets were a Aql (mV= 0.77 mag), a Boo (mV= -0.04 
mag), a Lyr (mV= 0.03 mag), /? Ori (mV= 0.12 mag) and a Cyg 
(mV*: 1.2 mag, variable star). 

Fig-Spresents the statistical distributions of the high angular 
resolution parameters. The seeing values (top left) span over the 
range [0.4 and 2.0 arcsec], with a half of the seeing better than 1.17 
arcsec, and 75% of the seeing better than 1.4 arcsec. The isoplanatic 
angle (top right) seems to follow a bimodal distribution with an 
average of 2.17 arcsec, and a median of 2.01 arcsec, respectively. 
The coherence time (bottom left) spans over a 1. to 5 ms range with 
an average and median values of 2.02 and 1.83 ms, respectively. 
For the coherence etendue (bottom right), the average and median 
values are 0.14, and 0.06 m 2 . ms. arc sec 2 . 

Fig- H] presents, night by night, the temporal evolution of the 
optical turbulence over the whole campaign, from ground level up 
to 30 km, along with detailed explanation given in the caption. 

Fig[6]shows the mean vertical profiles of turbulence and wind 
speed over the whole atmosphere. Almost all the optical turbulence 
is concentrated within the boundary layer, and the rest being scat¬ 
tered through the free atmosphere. The wind speed, shown on the 
top right panel of Fig(6] increases from 10 m.s~ l at ground level to 
reach 24 m.s~ l at 2 km. As discussed in Giordano et al. (2012), the 





















































Optical Turbulence Characterization at LAMOST Site 5 



Figure 5. Temporal evolution of the C^(h) during the 7 nights where SSS was observing. The upper part of each plot has a one-km vertical resolution and the 
lower part refers to a kind of zoom in the boundary layer (first km of the above plot). One can see that most of the time the first km (lower part of the upper 
frame) is reddish, corresponding to C 2 N ~ 10 15 m ~ \ Within the free atmosphere, 1 km above ground level, the optical turbulence is faint and scattered from 
1 to 30 km. When C 2 N profiles stops below 30 km, this means that the line of path of the star was far from zenith and, when straighten to the vertical, it was not 
possible to reach 30 km. See Habib(2006) and Vernin(2009) and references therein. 


SSS can retrieve four values of turbulence and wind speed within 
surface layer. Due to its low vertical resolution, the SSS technique 
is not able to distinguish the altitudes of the four layers. Typically, 
there is a wind gradient above a few hundred meters of Earth’s sur¬ 
face. Therefore, we arbitrarily sorted these layers with increasing 
wind speed from bottom to top of the boundary layer, which seems 
reasonable and already checked by Giordano et al. (2012). The bot¬ 
tom panel of Fig[6]shows the profiles of turbulence and wind speed 
modulus within the first four sub-layers, corresponding to a surface 
layer (0-1 km) slab. One can notice the steep wind speed gradient 
from 6 m.s~ l to 14 m.s~ l in the boundary layer. 

The relative error of the correlation function is defined by 
Vernin et al. (1983): 


AC 

C(0j 

where AC is the estimation error, N spec kie S is the amount of indepen¬ 
dant atmospheric speckles analysed, and C(0) is the scintillation 
variance which is created by all of the layers. As is well known, 
the scintillation index for a given turbulent layer with altitude h as 


N, 


speckles 


(6) 


expressed by Roddier (1981): 

A of = AC = l9.2A- ll6 h 5l6 A(C 2 N dh), (7) 


The SSS absolute error can be deduced with formula (6) and (7): 

A ll6 lr 516 


A(Qcft) = 0.052C(0)' 


iV speckles 


( 8 ) 


The scintillation index crj (the normalized intensity variance) 
is related to the C 2 profile (Roddier 1982), 


cr 2 = 19.2A- J/6 ^ h 5/6 C 2 (h)d/ij. (9) 

where I is the normalized intensity fluctuation. Here, the scintil¬ 
lation is measured on each point of the telescope pupil, which are 
k 0.7 cm spaced, much less than the first Fresnel zone which ~ 2cm 
for a 1 km altitude layer. One can assume that the pupil averaging 
is almost negligible. The autocorrelation of the scintillation C,{r) 
measured by the SSS on the pupil plane could deduce the scintil¬ 
lation index cr] = C»(0). Fig[7]is the statistical distribution of the 
scintillation index at LAMOST site in 2011; the mean cr] is 0.27, 
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-SSS- xl-40 - 2011 



Figure 6. The profiles of turbulence C„(h) and wind speed deduced by the 
SSS measurements at the LAMOST site in 2011. The top-left is the median 
profile of turbulence C„(h) with ±cri og (CN 2 ) error bars, and the top-right the 
median of the wind speed modulus profile \V(h)\ with ±<t error bars. The 
average profiles of turbulence (bottom-left) and speed modulus (bottom- 
right) within the surface layer (4 sub-layers) are also shown. 


Xinglong — Cn2 mean Xinglong — |V| mean 




Cn2 |V| (m.s 1 

Figure 9. The comparison of the average profiles of turbulence C%(h) and 
horizontal wind velocity \V(h)\ on 18 April 2011, obtained by the SSS, by 
WRF model, and by the NOAA, the Air Resources Laboratory of National 
Oceanic and Atmospheric Administration. The blue line by the NOAA 
shows the wind speed above Xinglong observatory at 21:00 UT on 18 April 
2011 . 


Statistical distribution of a. 2 with SSS in 2011 



Figure 7. Histogram and cumulative distribution (dashed line) of the scin¬ 
tillation index crj measured by the SSS at the LAMOST site in 2011. 


and the median crj is 0.25. That crj < 0.6 for 99% of the scintilla¬ 
tion index indicates the regime of weak fluctuations, and allows to 
use the small perturbation theory (Roddier 1981). 


ing is better than 1.4 arcsec. The range of the isoplanatic angle 0 0 
deduced from the WRF model is two times smaller than the SSS 
measurements, but the median and mean values are close (2.26 arc- 
sec and 2.25 arcsec with WRF model compared to 2.01 and 2.17 
with the SSS, and their cr is 0.78 arcsec). The distribution of the 
coherence time by WRF model is really different from the SSS 
measurements, and the median and mean values are 4.06 ms and 
3.90 ms, respectively, much greater than the SSS results. The co¬ 
herence etendue Go by WRF model has a log-normal shape similar 
as that by the SSS, and the median and mean Go are 0.18 and 0.35 
m 2 . ms. arc sec 2 , also close to the SSS values. 

Fig[9] presents the average vertical profiles of the C 2 (/r) as 
measured by the SSS and by WRF model, and the vertical pro¬ 
files of the average horizontal wind velocity \(h), by the NOAA, 
by the SSS, and by WRF model. It seems that the profiles by WRF 
model do not meet those by the SSS measurements. The C 2 (/i) pro¬ 
file by the SSS is distinctly larger than that by WFR model; the 
wind profile by the SSS shows higher velocity in lower atmosphere 
and lower around the tropopause, and looks quite different from 
that by WFR model. The vertical profile of wind velocity by WRF 
model and by the NOAA looks consistent, though there is about 1 
km offset in altitude, with a minimum at the ground level of 4 m.s~ l 
and a maximum around 12.5 km of 35 m.s~ l . 


3.3 WRF model results 

We used the WRF model to forecast the optical conditions at the 
same time as the DIMM and SSS run. The WRF model was con¬ 
figured to compute one forecast every 10 minutes in order to have 
a good temporal resolution. 

Fig[8] shows the statistical distributions of the four integrated 
atmospheric optical parameters. The seeing has a bimodal distri¬ 
bution with the first bump at 1.1 arcsec and the second at 1.6 arc¬ 
sec. The median and mean seeing values are 1.07 arcsec and 1.12 
arcsec, respectively, very close to the DIMM seeing values (1.01 
arcsec and 1.09 arcsec), the standard deviations cr is 0.28 arcsec. 
Moreover, the cumulative distribution shows that 75% of the see- 


4 DISCUSSION 

We have performed a campaign at the LAMOST site in 2011 to 
characterize the optical turbulence with the instruments of DIMM 
and SSS, and by the WRF model. The correlation analysis are per¬ 
formed with the results in order to compare the DIMM and SSS 
measurements, as well as to check for the WRF model previsions, 
which have been already validated by Giordano (2013). 

Comparing the seeing distributions as measured by the DIMM 
(Fig[3}. by the SSS (FigH and WRF model (Fig[H, one can notice 
that mean and median values are coherent, but the shape of the dis¬ 
tribution looks somewhat different. WRF seeing distribution shows 
a bimodal behavior whereas DIMM peaks at low seeing values and 
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Figure 8. Histogram (solid line) and cumulative distribution (dashed line) of atmospheric parameters forecasted with WRF model, from April to September 
2011 at the LAMOST site. 


SSS at large seeing values. If one assumes that free atmosphere op¬ 
tical turbulence is steady, as seen in Fig[5] and the boundary layer 
turbulence is more intermittent, one should expect a bimodal dis¬ 
tribution as already encountered by Aristidi et al. (2009) in Antarc¬ 
tica. But, we have no explanation to the different DIMM and SSS 
seeing distribution. 

Fig|T0]shows the comparisons of the seeing averaged every 10 
minutes by the three methods. Statistically, DIMM and SSS seeings 
agree within a standard deviation of cr=0.4 arcsec. WRF seems to 
predict larger seeing values than SSS and DIMM, with a regres¬ 
sion slope of 0.8 and 0.78. Even if the three regression slopes are 
close to one, meaning that there is no systematic bias, one can no¬ 
tice a large scatter. Indeed, the standard deviations of the linear 
fitting show cr=0.41 arcsec between the DIMM and SSS seeing, 
<x=0.44 arcsec for the DIMM and WRF seeing, and <x=0.45 arcsec 
for the SSS and WRF seeing. The correlation between the SSS and 
DIMM seeing seems to be a little better. For site characterization, a 
10 minute resolution seems too high, and a one day seems a more 
appropriate sampling. Fig[TT] further compares the seeing values 
averaged over each night. One can see that not only the regression 
slope is closer to one, but the dispersion a is almost two times bet¬ 
ter. If one assumes that SSS and DIMM are well characterizing the 
seeing, WRF seems a good tool for seeing characterization, as al¬ 


ready mentioned by Giordano et al. (2013) and discussed in section 

m 

The comparisons of the other turbulence parameters measured 
by the SSS and calculated by WRF model are shown in Fig[l2] for 
the isoplanatic angle do, the coherence time tq and the coherence 
etendue Go, and also the parameters for free atmosphere (FA): e FA , 
tfa and Gfa ■ We find that the agreement is not good; though for 
the seeing efa the deviation is small (cr=0.13 arcsec), for the To and 
Go, it is not the case. The isoplanatic angle is dominated by high 
altitude (FA) layers which are predicted within a small dynamic 
range by WRF. WRF optical turbulence relies on Trinquet-Vernin 
parametrization which is based upon vertical gradient of the poten¬ 
tial temperature and wind speed, but the vertical spacing is coarse. 
This means that FA dependent parameters such as Sfa and do are 
not subject to large variations. 

Tab[T|compares the optical turbulence parameters at the LAM¬ 
OST site with those at ORM, Mauna Kea, Armazones, and the 
Antarctica sites as Dome C and South Pole, (Vernin et al. 2011, 
Schock et al. 2009, Giordano et al. 2012, Marks et al. 1999). We 
find that the isoplanatic angle do at the LAMOST site is similar to 
the other sites, but the seeing and coherence time at the LAMOST 
site are poorer. The poor coherence etendue Go, at the LAMOST 
site, as well as South Pole, may be not good for adaptive optics. 
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SSS vs DIMM comparison in 2011 


DIMM vs WRF comparison in 2011 


SSS vs WRF comparison in‘ 2011 





Figure 10. The correlation maps of 10 minutes average seeing eo for SSS vs DIMM, DIMM vs WRF, and SSS vs WRF. The black line means complete 
correlation, and the red line is the linear fitting. 


SSS vs DIMM doily correlation in 20T;.t 


DIMM vs WiRF daily correlation in 2011 


vs WRF daily correlation in 201 1 





SSS seeing WRF seeing 

Figure 11. Same as figure flol but for the seeing values averaged over one night. 


WRF seeing 


This can be due to the high optical turbulence distribution in the 
boundary layer, associated to a steep wind shear, from 6 up to 14 
m.s~ l in the first four ground layers (See Fig[6] the bottom right 
panel). 

In Tab(2] we recall all of the site testing campaigns at Xinglong 
observatory, and we collect the seeing results every measurements 
in different place, in order to provide a reference to the application 
of Adaptive Optics(AO), the active optics and the high resolution 
spectrometer. Song et al. (1998) measured the seeing in 1994, using 
a DIMM inside the 2.16 m telescope dome, and the results could be 
affected by the dome seeing. Liu et al. (2003) and Yao et al.(2012) 
analyzed the FWMH of Polaris images obtained by a Schmidt tele¬ 
scope for a long period of time (1995- 2011). Although the results 
includes many influence factors, such as the dome seeing, image 
quality, and tracking accuracy of the telescope, however, the data 
collection for over a decade is useful to understand more about the 
seasonal seeing conditions. Liu et al. (2010) performed DIMM see¬ 
ing measurements on the top of LAMOST building in 2007, and 
obtained a better seeing with higher elevation than that at the 2.16m 
telescope. In this work, we selected an open flat roof for the mea¬ 
surements, and for the first time, we obtained whole night optical 
turbulence profiles; this work should provide more reliable results 
to characterize the site condition at Xinglong observatory. 


5 CONCLUSION 

A site testing campaign has been carried out in 2011 at the LAM¬ 
OST site in order to characterize the optical atmospheric turbu¬ 
lence. A DIMM and a Single Star Scidar have been set up during 
the campaign, and the WRF model is applied to calculate the site 
condition and compare with the observation results. For the first 
time, we are able to measure the turbulence profiles at Xinglong 
observatory. 

The results of the optical turbulence parameters - seeing, iso- 
planatic angle, coherence time, and coherence etendue, are pre¬ 
sented. The integrated seeing measured by the DIMM gives a me¬ 
dian value of 1.01 arcsec. There are 6011 individual profiles of 
optical turbulence C 2 (/z) and horizontal velocity |v(/i)| recorded by 
the Single Star Scidar, and the median SSS seeing is measured to 
be 1.17 arcsec, the median isoplanatic angle 2.01 arcsec, and the 
median coherence time 1.83 ms. The global parameter coherence 
etendue Go, which is usually employed to evaluate the performance 
of high angular resolution, gives a mean value of 0.14 and a median 
of 0.06 m 2 .ms.arcsec 2 . 

The optical turbulence parameters are calculated over the 
same period of the campaign with the WRF model, coupled with 
the Trinquet-Vernin model, based on the global weather data by Na¬ 
tional Center for Atmospheric Research (NCAR). The WRF model 
calculation gives a median seeing of 1.07 arcsec, a median isopla- 
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Figure 12. The correlation maps of 10 minutes average optical turbulence parameters measured with the SSS and forecasted with WRF model. The upper 
diagrams are the isoplanatic angle 60 , the coherence time to and the coherence etendue Go, and the below are the respectively parameters of free atmosphere 
(FA): the seeing spa, tpa and Gpa■ The black lines mean complete correlation, and the red lines indicate the linear fitting. 


Parameters 



Site 



ORM 

Mauna Kea 

Armazones 

Dome C 

South Pole 

Xing Long 

Total seeing £o (arcsec) 

0.80 

0.75 

0.64 

1.00 

1.60 

1.17 (SSS) 

Isoplanatic angle #o (arcsec) 

1.93 

2.69 

2.04 

6.90 

3.23 

2.01 (SSS) 

Coherence time to (ms) 

5.58 

5.10 

4.60 

3.40 

1.58 

1.83 (SSS) 

Coherence etendue Go (m 2 ms arcsec 2 ) 

0.38 

0.62 

0.49 

1.80 

0.07 

0.06 (SSS) 

Reference to data 

Vemin (2011) 

Schock (2009) 

Schock (2009) 

Giordano (2012) 

Marks (1999) 

This article 


Table 1. The global median values of optical turbulence parameters. 


Period 

Method 

Night 

Elevation(m) 

Spring(arcsec) 

Summer(arcsec) 

Autumn(arcsec) 

Winter(arcsec) 

Reference 

Dec. 1994 

DIMM 

4 

15 




1.4 

Song (1998) 

1995-2001 

FWMH 

All 

6 

3.5 

2.9 

3.7 

3.9 

Liu (2003) 

Oct.2003 

DIMM 

3 

15 




1.3 

Zenno (2004) 

Mai-.-Apr.2007 

DIMM 

12 

28 

1.1 




Liu (2010) 

Apr.-Sep.2011 

DIMM 

7 

4 

0.91 


1.09 


This article 


SSS 

7 

4 

1.32 


1.10 




WRF model 

8 

" 

1.12 


1.04 




Table 2. The seasonal comparison of the median seeing measured with multi-methods at different locations of Xinglong observatory. 
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natic angle of 2.26 arcsec, and a median coherence time of 4.06 ms. 
The coherence etendue Gg by the WRF model gives a mean value of 
0.35 and a median of 0.18 nr. ms. arc sec 2 . We find the WRF model 
can provides more accurate seeing forecast on one night average 
than ten minute average. 

If WRF previsions seems to depart from SS S/DIMM measure¬ 
ment, this could be due to the fact that the l.rlkm terrain model 
might be to coarse when compared to the steep slant of about 200m 
slope situated along a SW-NE direction ridge, at the SE border of 
XingLong observatory. 

The optical turbulence parameters obtained in the campaign 
are compared with other sites in the world, and the seeing mea¬ 
surements in history at Xinglong Station are reviewed. The Coher¬ 
ence time To (1.58 ms) and Coherence etendue Go (0.07 m 2 ms 
arcsec 2 ) at LAMOST site is poor, leading to the conclusion that no 
exceptional conditions should be expected for adaptive optics. The 
C 2 (/7) measurements presented in this work are evaluated over only 
a short period. More DIMM and SSS measurements are needed 
in the future for a comprehensive optical turbulence evaluation of 
LAMOST site. 
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